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The response of spin-2 small condensates to an external magnetic field B is studied. The param- 
eters of the interaction are considered as variable. The emphasis is placed on clarifying the modes 
of excitation caused by the quadratic Zeeman term. The theoretical method used is beyond the 
mean field theory. A set of eigenstates with the U{5) D 5*0(5) D 5*0(3) symmetry is introduced to 
facilitate the analysis. To obtain a quantitative evaluation on the response, the fidelity susceptibility 
and the B -dependent average populations of spin-components have been calculated. Mostly the 
particle number A*' = 30 is assumed. The effect with a larger or smaller A'' is also considered. It 
was found that the sensitivity of the response depends strongly both on the interaction and on the 
inherent symmetry. 
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I. INTRODUCTION 

The physical properties of Bose-Einstein condensates 
of atoms with nonzero spin can be tuned by external 
laser fields and magnetic fields T] . These condensates are 
much less affected by defects and impurities which exist 
extensively in usual condensed matter systems. There- 
fore, the condensates can be used to simulate a variety of 
ideal condensed matter systems (e.g., the system of ghi- 
ons) for academic studies, or to work as ideal magnetic 
materials in practical application (say, in quantum infor- 
mation and quantum computation). For the purpose of 
manipulation, how the condensates will react to an ex- 
ternal magnetic field is a crucial problem to be clarified. 

Usually, the particle number N in a condensate would 
> 10**. If could be greatly reduced, new physics might 
emerge. For an example, for spin-2 condensates, it was 
predicted based on the mean field theory (MFT) that the 
ground state would have three phases, namely, the fer- 
romagnetic phase (/-phase), polar phase (p-phase), and 
cyclic phase (c- phase) However, from a study based 

on a theory beyond the MFT, the ground state may have 
3, 4, or 5 phases depending on iV [Sj. The fine effect of 
N could only be seen if N is small (say, N < 100). When 
N is small, the properties of the condensates might de- 
pend on TV sensitively, just as the properties of the nuclei 
depend on the number of nucleons sensitively. On the 
other hand, recent progress in technique suggests that 
the condensates with N very small could be experimen- 
tally achieved @. Therefore, a study of the small con- 
densates is worthy because new knowledge additional to 
those from large condensates might be obtained. 
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This paper is dedicated to small spin-2 condensates, 
and is a generalization of two previous papers 0, Q on 
spin-2 and spin-1 condensates. The method adopted 
is beyond the MFT. The aim is to study the sensitiv- 
ity of the small condensates responding to the variation 
of an external magnetic field B. To this aim the fi- 
delity susceptibility [9l-[Tl| and the average population 
of spin-components have been calculated. To clarify 
the underlying physics, the set of eigenstates with the 
U{5) D 5*0(5) D 5*0(3) symmetry firstly proposed in 
the refs. [1, [l^ is introduced. Based on the set, elemen- 
tary excitation modes caused by the quadratic Zeeman 
term have been found. These modes together with the 
associated energy gap are decisive to the response of the 
spin-structures against the appearance of B. It is well 
known that the structures of the ground states depend 
on the interaction. On the other hand, the ways of ex- 
citation caused by the field is found to be strongly con- 
strained by the inherent symmetry. Thus both the inter- 
action and symmetry together determine the response of 
the condensates to the field as shown below. 

Since the study is not dedicated to a specific spin-2 
system, the parameters of the spin-dependent interaction 
are considered as variable. Mostly TV = 30 is given, the 
effect of a larger or smaller N is also discussed. 



II. HAMILTONIAN 



The atoms are assumed to be confined by an isotropic 
and parabolic potential imw^r^. The temperature T is 
assumed to be very low. Since the energy of spatial exci- 
tation of individual atom is ~ huj, when T << hu/kB it 
is reasonable to assume that no atoms would be spatially 
excited and all of them would fall into a common spatial- 
state 4>{r). This assumption is called the single mode 
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approximation (SMA), which has been quite frequently 
used in the hteratures. Furthermore, a magnetic field B 
lying along the Z-axis is set. When an integration is car- 
ried out over the spatial degrees of freedom, we arrive at 
a model Hamiltonian depending only on the spin-degrees 
of freedom as[^ Il2l | 



(1) 



where an overall constant has been neglected, Vij = 
ds^'Ps ' where S is the combined spin of the particles 
i and j, Vg is the projector of the S'-channel, gs is the 
strength of interaction related directly to the s-wave scat- 
tering length of the S'-channel, n = J \cj){r)\'^dr, and the 
last two terms are for the linear and quadratic Zeenian 
energies, respectively. Note that V!^^ = I — Vq —1^2- We 
further introduce l3 and 9 so that l3 cos 9 — {go — gi) / Qi, 
Ps\n9 ^ (52-54)/ff4- Then we define H'^^^ = jj^H^od- 
Furthermore, since M (the Z-component of the total 
spin) remains to be conserved under the field B, the lin- 
ear Zeeman term provides only a constant and therefore 
can be neglected. When all the constants involved are 
removed, -ff^iod '^^^ rewritten as 



^mod-'E(cOS^P^^ 



sm9r^) + q'Y.if.zf, (2) 



where 9 is from to 27r, q' = q/{(3g4n). When the Hamil- 
tonian is transformed from Eq. (IT|) to ([2]), the set of 
eigenenergies will be multiplied by a common constant 
and will be further shifted as a whole. But the eigen- 
states will remain unchanged. Since only two parameters 
9 and q' are contained in H'^^^, related analysis is easier 
to perform. 

When q' — 0, it has been proved that H^^^ can be 
rewritten as a sum of four Casimir operators of a chain 
of groups C/(5) D 50(5) D 50(3) (3. Therefore the 
problem can be solved analytically based on the group 
algebra [5l. [isl. . Three quantum numbers, namely, F 
(total spin), M (its Z-component), and v (seniority, N—v 
must be an even integer, and (N — v)/2 is the number of 
singlet pairs) are introduced to classify the states. The 
eigenenergies related to these quantum numbers are 

^ 10 sin — 7 cos , ^ am9 , ^ , 
E.F^ v{v + i)--^F{F + l), (3) 

where an irrelevant overall constant depending on N has 
been further neglected. The associated eigenstates are 
denoted as '^vFM- Note that, for specifying a symmetry- 
adapted eigenstate, not all the (w, F) pairs are allowed. 
Some of them are prohibited by inherent symmetry. The 
details are referred to (13.]. Eq. ([3|) leads to a notable 
feature, namely, when s\a9 = ^ cos 6*, all the symmetry- 
adapted '^vFM with the same F are degenerate. Whereas 
when sin = 0, all the symmetry-adapted '^vFM with the 
same v are degenerate. If the group of degenerate states 
happens to be the ground states, the low-temperature 
behavior of the condensate would be seriously affected. 



In general, a slight increase of N (say from being even 
to odd) will affect the structure of the ground state . 
However, when N is large, the effect is weak. For the con- 
venience of discussion, it is assumed that is a multiple 
of 2 and 3 in the follows (other cases of N can be similarly 
discussed). When 9 and M are fixed [M > is assumed), 
from Eq. Q and from the symmetry-adaptability, the 
pair {vg,Fg) for the ground state (the lowest state with 
the given M) can be determined. When M = 0, based on 
the method of group algebra, there are three regions of 
9 5]. When 9 is from to 9fp = arctan[— iqIJ^_^2) ] 
(region V), Vg — N and Fg = 2N, and accordingly the 
ground state is in the /-phase. In this phase all the spins 
are nearly aligned along a common direction to form a 
total spin, which is vertical to the Z-axis (if M — 0) and 
has an arbitrary azimuthal angle. When 9 is from 9fp 
to 9pc = arctan(^) = 214.99° (region II), Vg = and 
Fg = 0, and accordingly in the p-phase, in which the 
spins are two-by-two to form the singlet pairs (If N is 
odd, the ground state would have Vg = 1 and Fg = 2, 
and accordingly would be composed of a group of singlet 
pairs together with a single particle. This case is not dis- 
cussed in this paper ). When 9 is from 9pc to 27r (region 
III), Vg = N and Fg — 0, and accordingly in the c-phase, 
in which the spins are essentially three-by-three to form 
the triplexes (If is not a multiple of 3, additional parti- 
cle(s) would be added) The division into three regions 
based on the group algebra coincides with that from the 
MFT when A^ is sufficiently large[2j, |4|. Nonetheless, the 
details of spin correlations can be understood from the 
former but can not from the MFT. 

When q' 0, there is no analytical solution. F 
is no more conserved, but M remains to be a good 
quantum number. To obtain numerical solution, the 
Fock-states |a) = | A^a": ^i": ^o": ^-i- ^-2) are used as 
basis functions for the diagonaiization of Ft^^^, where 
A^^ is the number of particles in the /x-spin-component, 
= N and f^N^ = M, \a) as a whole form a 
complete set for symmetrized spin-states. Thereby exact 

eigenstates V'^m '^^ ^'mod *^an be obtained, where z is a 
serial number of the series of states with the given M 
(i = 1 denotes the ground state). Various information 

will be extracted from V'o//- '^^se with M = is 

firstly considered. 



III. THE FIDELITY SUSCEPTIBILITY OF THE 
GROUND STATE AND THE MODE OF 
EXCITATION 



Firstly, we would like to study the sensitivity of the 
ground state V'im against the variation of the field. For 
this purpose the fidelity susceptibility [§-[ll| 



XM(<z') = lim^(l-|(^?;rWM)|), 



(4) 
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FIG. 1: (Color online.) ln[xo(0)] (solid) and HEg^p) (dash) 
against 9. AI — is assumed. The two vertical dotted lines lo- 
cated at 9fp and 6pc, respectively, are introduced for separat- 
ing the regions I, II, and III. In (a) iV = 30 and Ofp = 140.48° , 
while in (b) iV = 6 and Ofp = 122.41°. In both (a) and (b), 
dpc = 214.99°. 



has been calculated. A larger xmW) implies that the de- 
viation between V'im'^ and V'im larger, thus the state 
responds more sensitively to q' . Therefore this quantity 
measures quantitatively the sensitivity. From the eigen- 
states V'l and V'l o (both can be obtained via the diago- 
nalization), ln[xo(0)] of the ground state as a function of 
9 is shown by the solid curves in Fig. [T] The behaviors 
of Xo(0) in the three regions of 9, which are separated 
by the dotted vertical lines, are different. In general, the 
p-phase in region II is much more sensitive to the ap- 
pearance of the magnetic field than the /- and c-phases. 
From Eq. ([3]) , we know that the ground state is degener- 
ate when 9 = 9pc or 0. In these cases, the susceptibility 
can not be well defined. Therefore, in Fig.[Tl the two nar- 
row domains {9pc — S, 9pc + 6) and {—6, 6) with S — 7r/180 
are actually not included in the calculation. Nonetheless, 
when 9 is close to these narrow domains, the sensitivity 
is very high. 

In order to understand the behavior shown in Fig. [TJ 
we have to study the mode of excitation caused by the 
quadratic Zeeman term. The rule of selection governing 



the matrix elements 



Qv'F',vF — {'^v'F'M\y^Xfiz)'^\^vFM) 



(5) 



is crucial. One can prove that Q^'p, is nonzero only 



if 



= or ±2 and \F' - F\ < 2. Additionally, if 



M = 0, (—1)^ = (^1)^ is required. Of course, both the 
pairs {v',F') and {v,F) should be symmetry-adapted. 
Let the excitation caused by Q^lp, yp be called the ele- 
mentary excitation. For the ground state ^vgFgO, due to 
the symmetry constraint [isj. one can prove that the non- 
diagonal Q^ipi y p is nonzero only if {v',F') = {vg,Fg), 
where {vg,Fg) = {N, 2N - 2), (2, 2), and {N - 2, 2), re- 
spectively, for the /-, p-, and c-phases (refer to Tab. U). 
It implies that each phase has only a unique mode of el- 
ementary excitation. This is a distinguished feature of 
M = states. 

When the ground state is not close to being degener- 
ate, the first-order perturbation theory can be used to 
calculate the susceptibility. When q' — and e is very 
small, due to having only a unique elementary mode, the 
first order approximation of ipf q can be written as 



i'U = li'^vgFg,o+e- 



(6) 



where i?gap — p — E^^ p^, 7 is simply a constant for 
the normalization and can be easily obtained. Inserting 
Eq. ^ into (g]), we have 



Xo(0) « (- 



E, 



(7) 



From Eq. we know that, for the /-phase Sgap = 
En,2N-2 - Em,2N = sin 9, for the p-phase E^gap = 

E2.2 — Eo_o = sin9 — cos 9, and for the c-phase E^gap — 

77 IT 7 cos 6^ — 10 sin j\7- I r,\ 3 sin ^ jji 
rjN-2.2 — J^N.n ~ (4iV +2) = — . iigap 



are also plotted in Fig. [1] by the dash curves. In general, 
Eq. ([71) gives a very good approximation except that 9 
is close to 9pc or (27r). When 9 varies within a re- 
gion, p p does not depend on 9. Therefore the 

6'-dependence is due to the factor (l/i?gap)^ and can be 
shown by comparing the solid and dash curves. Inciden- 
tally, since '^vgpg might not be the first excited state, 
the gap defined here might not be the energy difference 
between the ground and the first excited states. Further- 
more, when q' is not very weak, one has to go beyond the 
first order perturbation theory, and therefore, in addition 
to the elementary excitation, higher order excitations will 
be included. Accordingly, various ^t,',_F',o components 
will be contained in ■0^ q . 

The features of the ground states with M = are 
summarized in Tab. U] Some details in the table may be 
changed if N is not a multiple of 2 and 3 (Say, if N is 



-7(Af+4)i 



< vr 



odd, 9fp is changed to a.Tctaii[- ^g^j^_^^ 

Note that the elementary excitation of the /-phase in 
region I, namely, from {vg ^ N,Fg = 2N) to {N, 2N-2), 
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TABLE I: The three phases of the ground states with M = and the associated three connected regions of 6. N is assumed 
to be a multiple of 2 and 3. {vg,Fg) are the quantum numbers of the ground state, {vg,Fg) specifies the corresponding 
unique mode of elementary excitation, and iSgap is the associated energy of excitation. f v f '^^ matrix element of the 

elementary excitation caused by the quadratic Zeeman term This matrix element is given at A'' = 30 and 6 (inside 

the parentheses). 



region 



I (ferromagnetic) 



II (polar) 



&pc — arctan(TjZj) > vr 

0,0 

2,2 

sin 9 — cos 6 
28.98 (7.27) 



III (cyclic) 



right border 

Vg,Fg 
^g,Fg 
E 

o° - 



arctan 



-7(iV+3) i 
10(]V-2) J 



< IT 



N, 2N 
N,2N -2 

6.68 (2.93) 



Ocf = 2tt 
N,0 
N -2,2 

7 COB 9 -10 sin 9 
70 

6.93 (3.10) 



(4iV + 2) 



3 sin £ 
7 



is essentially a change of the total spin. The associated 
energy gap is nearly proportional to TV. Due to the en- 
largement of the gap by the factor N, the /-phase has a 
very low susceptibility in general. However, when — >■ 0, 
the gap tends to zero resulting in a great increase in the 
susceptibility. This explains the sharp peak appearing in 
the left end of Fig. [1^ or [TJj. Nonetheless, since the gap 
will increase with N, the sharp peak will be more and 
more suppressed when N is larger and larger. 

The elementary excitation of the p-phase in region II, 
namely, from (0, 0) to (2,2), is a transformation of a sin- 
glet pair into a S* = 2 pair. The associated (522. oo 

roughly proportional to N [l6|. Even when N is not 
large, it is still much larger than the p p of I and 

III as shown in Tab. H] Furthermore, the gap is small 
and does not depend on N . Therefore, the p-phase has a 
very high susceptibility, and will become higher when N 
is larger. The minimum of the gap appears at the right 
border of II. Accordingly, there is a peak [l7| . 

The elementary excitation of the c-phase in III, namely, 
from {N, 0) to {N — 2, 2) is essentially a transformation 
of a triplex to a singlet pair plus an extra particle. The 
associated gap is also enlarged by N. It will be in general 
large, therefore the c-phase has also a very low suscepti- 
bility. However, the gap would decrease very fast if 9 is 
close to 9pc. Accordingly, there is a sharp peak at the left 
border of HI as shown in Fig. [1] Up to now the feature 
in Fig. [T] has been explained. 

Comparing Fig. [1^ with [TJd we know that the quali- 
tative feature of the susceptibility does not change with 
N. However, since 9fp depends on N, the border be- 
tween I and II will move a little to the right when N be- 
comes larger and will tend to arctan(— 7/10) = 145.01° 
when N ^ oo (this value coincides with that from the 
MFT). Furthermore, the increase of N will cause a re- 
markable increase of Xo(0) in II (due to the fact that 
Q2 2 is roughly proportional to N) and will cause a 
small decrease of Xo(0) in I and HI (Although Q° f, „ 
increases with N , however the gap increase also. The 
two effects compete with each other, and finally the lat- 
ter over takes the former). 



IV. POPULATIONS OF SPIN COMPONENTS 
UNDER A MAGNETIC FIELD 

When a magnetic field is applied, since the quadratic 
Zeeman term imposes additional energy to the particles 
with /i 7^ 0, the ground state will prefer to have the num- 
ber of /J, = particles, iVo, larger so as to reduce the 
energy. Thus, in addition to the susceptibility, the study 
of the variation of Nq is also a way to understand the ef- 
fect of the field. Note that the experimental measurement 
of the susceptibility might be difficult. On the contrary, 
A^o can be easily measured, say, via the Stern-Gerlach 
technique. 

Let the eigenstates of i?4iod be expanded in terms of 

the Fock-states as = c^**^!*^)- One can extract 
a particle (say, particle 1) from a Fock-state as 

/ Kla 

where r\i^ is the spin-state of a particle in the component 
I', I • • • , A^^ - 1, • • • ) is a Fock-state of the [N - l)-body 
system, in which the number of particles in the compo- 
nent V decreases by one. Therefore, the particle can also 
be extracted from the eigenstate as 

^ Y.^Ami'^\ (9) 

V 

a 

Thereby the probabihty of a particle in u is just 

p^'^-^ ^ (V^f^Vf'")- (11) 

They fulfih Pu^^ = 1- Pi;^' is called the 1-body prob- 
ability, and NPl'^ = is just the average population 
of the v component of the z-th state. 

The one-body probabilities Pq '^ of the ground states 
with M = and q' given at a number of values are plotted 
in Fig. [21 When q' = (solid curve), the ground state 
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has three choices of phases as mentioned. The eigen-spin- 
state of the /-phase is 



2N,0- 



(12) 



where, the special way of spin couphng (i.e., the com- 
bined spin of an arbitrary group of j particles is 2j) 
assures that the spin-state written in Eq. (fT2)) is nor- 
malized and symmetrized. It is straight forward to ex- 
tract a particle from the form of Eq. ([T^. thereby the 
1-body probability of the /-phase can be obtained as 
Pj^'" = (C2^.'°2Ar-2 where the Clebsch-Gordan co- 

efficient has been introduced. When ly — 0, Pq'^ = 

(4N-i)^(4N-2K'iN-3) ■ ^his value is equal to 0.381 if 
TV = 30 as shown in Fig. H and will tend to 3/8 = 0.375 
if TV ^ oo. 

Both the p- and c-phases have F — {ii M — 0), and 
therefore no special orientation is preferred. Accordingly, 
as shown in = 1/5 for all v disregarding how N 

is. When the field is applied, both the /- and c-phases are 
inert to the increase of g', while the p-phase is extremely 
sensitive. These coincide with the previous findings in 
the susceptibility. In particular, when is close to the 
right border of II, a great increase of TVq appears even if 
the field is still weak. In I and III the increase of Nq is 
very slight except when is close to {2tt) and 6'pc, where 
the ground state is nearly degenerate as mentioned. The 
great difference in sensitivity with respect to q' holds 
disregarding how N is. Since the interactions of some 
realistic spin-2 atoms (say, ^'^Rb, ^^Na, and ®^Rb 0) 
might have 9 « 0pc, the high sensitivity in this particular 
domain would be very helpful to the exploration of the 
interactions. 
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FIG. 2: (Color online.) Pr 
M = 0. The solid, dash 
have q' — 0, 0.01, 0.02, and 0.03, respectively. 



Q against 9 with = 30 and 
dash-dot and dash-dot-dot curves 



Incidentally, if N increases from 30 to oo, the general 
feature of Fig. [2] would remain unchanged. However, the 
border between the regions I and II would shift from 
140.48° to 145.71°, the sohd line in I (for q' = 0) would be 



lower a little from 0.381 to 0.375, the parabolic-like curves 
in II for the p-phase would lie remarkably higher (say, for 
Af = 0, 61 = TT, and q' = 0.03, '° would be 0.282, 0.485, 
and 0.651 if = 6, 30, and 60, respectively). 
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(b) /V= 6, 
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FIG. 3: (Color online.) Pq against q' with = 30 (a) and 
6 (b), and M = 0. The solid, dash, and dash-dot curves have 
9 — tt/2, tt, and 3n/2, and represent the ferromagnetic (/), 
polar (p), and cyclic (c) phases, respectively. 



It was found that, when q' —5- oo, 

lim(0, 0, A^, 0, 0|f/'?o)^ = 1 disregarding how 9 is. It 
implies that all the ground states will tend to the same 
Fock-state |0,0,Af, 0, 0) as the strong field limit, where 
all the spins have ^ = 0. Accordingly, — s- 1 (if 
^ = 0) or — 7> (if ^ 7^ 0). The process of going to the 
limit is shown in Fig. [31 Where 9 is given at three values 
7r/2, TT, and 37r/2 to represent the /-, p-, and c-phases, 
respectively. One can see that the dash curve for the 
p-phase has a very steep take off at the left end. It 
implies that the p-phase needs only a relatively much 
weaker field to push it to its limit. This coincides with 
the finding from Fig. [2] When N is larger, the take off 
is steeper (e.g., when q' increases from to 0.25, P^'" 
would increase from 0.200 to 0.783, 0.877, and 0.905, 
respectively, if TV = 6, 30, and 60). It implies that, the 
p-phase of a larger condensate would be more sensitive 
to the field as found before. 

On the other hand, the take off of the dash-dot curve 
(c-phase) is much milder, and the solid curve (/-phase) is 
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the mildest, and they would be even milder if N is larger. 
For examples, when q' = 0.25, the dash-dot curves would 
be 0.300, 0.237, and 0.230, and the solid curves would be 
0.478, 0.393, and 0.384, if TV = 6, 30, and 60, respectively. 
It implies that, the /- and c-phases of a larger condensate 
would be more inert to the field than a smaller one. 



V. THE CASE WITH NONZERO 
MAGNETIZATION 

In the previous sections essentially the case with M — 
is considered. When q' = and M > 0, based on 
Eq. ([3]) together with the constraints 2v > F > M and 
(— 1)^ = (—1)", the dependence of the quantum numbers 
of the ground state {vg,Fg) on 9 are listed in Tab.HIl The 
effect of M is embodied by an integer vq, which is the 
smallest even (odd) integer > M/2 if N is even (odd), 
as shown in the table. Instead of three, there are four 
phases. However, when M = Ak (fc = 0, 1, 2, • • • ) and N 
is even, or M = 4fc + 2 and N is odd, we have 2vq = M. 
In this case, the two phases polar-A and -B are combined, 
and the four phases reduce to three as before. 

Note that, when M is < 5, (vo,M) and (iV, M) 
might be prohibited by symmetry [l^. Once {vo,M) 
or (N, M) is prohibited, other more advantageous 
symmetry-adapted (i;, F) pair will replace it (Say, in IIb 
where a smaller v and a smaller F will lead to a lower 
energy, if {vo,M) is prohibited, it would be replaced by 
{vo, M + 1) and/or (wq -|- 2, M). In III where a larger v 
and a smaller F will lead to a lower energy, if (A^, M) is 
prohibited, it would be replaced by (iV — 2,M) and/or 
{N, M + 1)). In these cases the appearance of more than 
four phases is possible. 

Just for an example, we choose N =even and M — 7 
(this choice is rather arbitrary). Then, we have vq = 4. 
Note that, for the region IIb, {vg.Fg) = (4,7) is prohib- 
ited. Therefore, it would be replaced by (4,8) and/or 
(6,7). Note that, from Eq. ([3]), we would have £'4^8 — 
Ee^r a e = arctan(^) = 208.16°. Accordingly,' the 
IIb splits into two. When n < 6 < arctan(:^), the 
symmetry-adapted ground state will have Vg = 4:,Fg = 8. 
Thus this part becomes an extension of IIa , and therefore 
the right border of IIa extends from tt to arctan(r^). 
When arctan(^) < 9 < arctan(j^) = 214.99°, Vg = 
Q,Fg = 7. Thus the region IIb with the quantum num- 
bers (6, 7) becomes very narrow. 

The fidelity susceptibility against 9 is shown in Fig. |4j 
where the division into four regions is clear (including the 
very narrow region IIb). In region I for the /-phase, the 
eigenstate '^n,2N,m = ({[(w)4?7]6?/}8 ■ • • )2Af,A/, where 
the spins are nearly lying along a common direction to 
form a total spin, which has an arbitrary azimuthal an- 
gle but not lying in the X-Y plane. In IIa and IIb, 
the ground states are still dominated by the singlet pairs 
together with unpaired particles. In IIa, the subsys- 
tem of the unpaired particles is in the /-phase, i.e., they 
form a spin-state as ({[(????)4f?]6f?}8 • ■ • )2vo,m- I'^ ^b, the 
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FIG. 4: (Color online.) In(x7(0)) against with iV = 30 
(solid) and 12 (dash). M = 7 is assumed. The domain of 6 
are divided into four regions separated by the dotted vertical 
lines (The dotted line at the left is simply the average of the 
two cases with A'' = 30 and 12). The quantum numbers of 
the ground states (vg,Fg) in each region are marked. 

spins of the unpaired particles are slightly diffused from a 
common direction, because their total spin M is slightly 
smaller than 2wo- In III, the c-phase is dominated by the 
triplexes together with a few unpaired particles. One can 
see that the susceptibilities of the pA- and pB-phases are 
in general much higher, and would become even higher if 
N is larger. This is similar to the findings from Fig.[TJ In 
particular, when 9 is close to the borders of the regions 
(except the border between I and IIa), the susceptibility 
is exceptionally high. 
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FIG. 5: (Color online.) P^'^ against 9 with M = 7 and 
N = 30. The solid, dash, dash-dot and dash-dot-dot curves 
have q' = 0, 0.02, 0.04, and 0.06, respectively. 

When q' is given at four values the populations of /i = 
components against 9 are shown in Fig. [5j Similar to 
Fig. [21 the pA- and pB- phases respond to the increase 
of q' more sensitively. The set of curves in IIa and IIb 
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TABLE II: The four phases of the ground states with M 0. Accordingly, the scope of 9 is divided into four connected regions, 
the right border of each region is given in the last row. vq is the smallest even (odd) integer > Af/2 if A*' is even (odd), and 
{vg, Fg) are the quantum numbers of the ground states. 



region 


I (/-phase) 


IIa (pA-phase) 


IIb (p-B-phase) 


III (c-phase) 


iV!„Fg) 

right border 


{N, 2N) 

arctani io[jv(]v-2)-t,o{"o-2)l ' 


(uo, 2«o) 

TT 


{vo,M) 

arctan(^) > vr 


(TV, M) 
27r 



will become even higher when N is larger. In particular, 
when 6 is close to {2tt) , and close to the narrow domain 
IIb , the sensitivity is very high. Obviously, the sensitivity 
emerging in Fig. [5] is closely related to the susceptibility 
shown in Fig. 21 

It was found that the number of nonzero non-diagonal 
matrix elements 0*^;=, would be in general more 

than one when M ^ 0. Therefore, there are more than 
one elementary modes of excitation (i.e., Eq. 1^ will con- 
tain more than two terms). The actual mode is a mix- 
ing of them. It was found that, when q' is very small, 
the ground state ipiM region I (/-phase), similar to 
Eq. ([6]), is essentially ^ N,2N,M + b'i' n,2N-2,m where 6 is a 
small constant. For a numerical example, when N — 30, 
M = 7,q' = 0.005 and 9 = 7r/2, b = -0.0019. This value 
will become larger when 6 is close to zero (say, b = 0.1072 
if = 7r/180). In general, |&| will decrease a little when 
N becomes larger. In the region IIa (pA-phase) and for 

M = 7 as an example, ^ ^ ^'4,8,7 + J2f^f^s,f,t^ 
where F is ranged from 7 to 10 and bp are in the order of 
1/10 or smaller. When 9 is close to the left border of IIa, 
only |6io| is relatively larger (meanwhile sin 9 is positive 
and the components with a smaller F will have a higher 
energy and therefore can be neglected, refer to Eq. ([3])). 
On the contrary, when 9 is close to the right border, only 
by I is relatively larger (meanwhile sin 9 is negative). In 
IIb (pB-phase) ^pfy « "Sajj + 6^*6, 8,7 + bp'i'sjj- In 

III (c-phase) V?'/ « *Af,7j + &7*Ar-2,7,7 + ^a^w.sj- The 
related coefficients in the two expansion are all small. 
When 9 is close to the left border of III, the gap asso- 
ciated with the change of v from TV to — 2 is very 
small, and therefore \b-y\ >> \bs\. On the contrary, when 
9 is close to the right border 27r, the gap associated with 
the change of F from 7 to 8 is very small, therefore 
|6^| « \bs\. 

When q' is large, it was found that limq'_j.oo (0, M,N — 
M, 0, 0|i/;^^)^ = 1 disregarding how 9 is, i.e., all the 
phases tend to the same Fock-state |0, M, N — M, 0, 0). 
Accordingly, P^'^'^ ^ M/N (if ^ = 1), (TV - M)/N (if 

fi — 0), and (otherwise). As examples, P/^'^and P^'* 
against q' are shown in Fig. [S] and Fig. [3 respectively. In 
both figures the solid curves ( /-phase) are inert to q' , the 
dash-dot curves (c-phase) respond to q' mildly. However, 
the dash curves (pA-phase in Fig. [5] and p-phase in Fig. [7]) 
are more sensitive to q' . In particular, at some critical 
points, abrupt changes are found in the dash curves im- 
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FIG. 6: (Color online.) against q with M = 7 and N = 
30. The solid, dash, and dash-dot curves have Q — it/3, tt, 
and 5n/3 to represent the /-, pA-, and c-phases, respectively. 
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FIG. 7: (Color online.) Pj'* against q' with M = 4 and 
TV = 30. The implications of the curves are the same as in 
Fig. [6] They are chosen to represent the /-, p-, and c-phases. 

plying a sudden transition in spin-structure (this abrupt 
change is not found in the previous cases of M — Q). In 
order to understand the transition better, the associated 
ground state is expanded as "01 m ~ 'Ylv F^fF^^FM- 
Then, we define Wy = J2f(^'vf)'^ which is the weight 
of having v unpaired particles. For the dash curve in 
Fig. [HI we found that W4 = 1 at = 0. When q' in- 
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creases, W4 decreases and We increases gradually. When 
q' — ^ q'^_^ = 0.098 (the critical point of the first tran- 
sition), W4 « 0.55 and Wq w 0.39. However, when q' 
crosses the critical point, W4 sudden becomes nearly zero 
while We jumps up to 0.90. Thus the transition is char- 
acterized by a great decrease of Wa together with a great 
increase of Wq. In Fig. [71 the dash curve have W2 = 1 at 
q' = 0. The first and unique critical point q'^_^ = 0.050. 
Similarly, the associated transition is characterized by a 
great decrease of W2 together with a great increase of 
W4. When q' is larger than the range of Figs. [S] and [71 no 
further critical points are found. Instead, all the curves 
tend smoothly to their limit M/N. 

VI. FINAL REMARKS 

The effect of a magnetic field on the spin-structures 
of the ground states of small spin-2 condensates is stud- 
ied. The elementary modes of excitation caused by the 
quadratic Zeeman term are found. The fidelity suscep- 
tibility and the populations of spin-components are cal- 
culated. Mostly the case = 30 is studied in detail. 
The effect caused by increasing and decreasing TV is also 
discussed. The following points are mentioned. 

(i) The set of eigenstates "^/vFM with the C/(5) D 
50(5) D SO{3) symmetry are introduced to help the 
analysis. When the magnetic field is zero {q' — 0) and 
M = 0, the pair of quantum numbers of the ground states 
(vg,Fg) have three choices associated with the three 
phases (/-, p-, and c-phases).0 It turns out that, for 
each phase, among all the non-diagonal matrix elements 
of the quadratic Zeeman terms, only the one « 
is nonzero. It implies each phase has its unique mode of 
elementary excitation. Thus the response of the system 
to a weak field is clear. In particular, the variation of the 
fidelity susceptibility against 9 can be understood in an 



analytical way. 

(ii) For the case Af 7^ 0, the ground state in general 
has four phases, namely, the /-, pA-, pB-, and c-phases, 
each has in general more than one elementary excitation 
modes. However, when 9 locates at some particular do- 
mains, one elementary mode will be dominant. 

(iii) The calculation of the fidelity susceptibility 
against 9 confirms the existence of the phases. The p-, 
pA- and pB-phases have in general a much higher suscep- 
tibility. In particular, when 9 is close to the borders sep- 
arating the phases, the susceptibility may be very high. 

(iv) In addition to the fidelity susceptibility, the vari- 
ation of the populations against the field has also been 
calculated. Both quantities provide similar message. In 
particular, the high sensitivity in the neighborhood of the 
borders is found again in the population. Since it has 
been suggested that a number of realistic species would 
have their parameters of interaction close to the borders, 
since the population is easier to measure, the high 
sensitivity is useful to the accurate determination of the 
parameters of interactions. 

(v) when g' — > 00, all the ground states tend to the 
limit |0, M, N — M, 0, 0) disregarding how 9 is. The pro- 
cess of tending to the limit is in general smooth. However, 
when M and q' is not very large (refer to Fig. [6]) , 
abrupt changes against q' are found. It implies that the 
field will cause not only a smooth variation of structure 
but also a phase transition at some critical points. 
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S is the operator of symmetrization, and the spins are 
two-by-two coupled to zero. ^'220 is exactly proportional 
to S{titi)2()[{'I'I'''I)(>]'^^'^^'^ ■ Since they have a rather simple 
form, they can be easily expanded in terms of the Fock- 
states. Thereby the matrix element for the elementary 
excitation, Q22,00) can have an analytical from. So, we 
can obtain numerical result even when N is larger (say, 
N = 500), and we find that (522,00 is roughly proportional 



to 

[17] In fact, when 9 is close to Opc., not only the gap decreases, 
a group of states having different v but the same F = 
Fg fall down and the ground state becomes to be nearly 
degenerate. This is also an important factor to raise the 
susceptibility. 



